Resonant enhanced diffusion in time dependent flow 
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Abstract 

Explicit examples of scalar enhanced diffusion due to resonances between differ- 
ent transport mechanisms are presented. Their signature is provided by the sharp 
and narrow peaks observed in the effective diffusivity coefficients and, in the absence 
of molecular diffusion, by anomalous transport. For the time-dependent flow con- 
sidered here, resonances arise between their oscillations in time and either molecular 
diffusion or a mean flow. The effective diffusivities are calculated using multiscale 
techniques. 

PACS number(s): 47.27Qb ; 47.27-i. 

1 Introduction 

Passive scalar transport in a given velocity field is an issue of both theoretical and ap- 
plicative relevance |l], |2|. The combination of molecular and advective effects can lead 
to rather subtle behaviors, even in the case of simple laminar flow. A rather complete 
theory has been developed for the dynamics at long times and large scales (macrodynam- 
ics), using multiscale techniques ||. For incompressible velocity fields the macrodynamics 
in the presence of scale separation is governed by an effective equation which is always 



diffusive and transport is always enhanced [|J. A sufficient condition ensuring the pres- 
ence of scale separation for three-dimensional incompressible flow is a vector potential 
finite variance, both for static and time-dependent flow || |6|] (with similar conditions 
for other space dimensionalities). The small-scale velocity field properties emerge via the 
effective diffusivity second-order tensor. Its calculation is reduced to the solution of one 
auxiliary partial differential equation that can be efficiently solved numerically @, || and 
variational principles have also been derived || . The effective diffusivities for generic flow 
depend on all the mechanisms of transport, e.g. on the advecting velocity field, on the 
molecular noise, on the presence and the strength of possible mean flow |7|, [TIJ. Each 
process is characterized by its typical time-scales and, if the latter are properly tuned, 
their mutual interference can produce resonance effects. These lead to strong enhance- 
ment of transport, reflected in sharp and narrow peaks in the effective diffusivity. Our aim 
here is to present explicit examples of such resonance effects in simple time-dependent 
incompressible velocity fields. In Section 2 the source of resonance is the presence of a 
large-scale velocity field. For random parallel flow (Sec. 2.1), explicit analytic expressions 
for the effective diffusivities are obtained and they indeed display strong peaks for specific 
strengths of the large-scale velocity. In Section 2.2 the robustness of resonance effects 
to variations of the scale separation ratio between the large- and the small-scale velocity 
field components is investigated. Numerical simulations show that the effects persist even 
for moderate scale separations. Section 3 is devoted to a time-periodic 2D flow widely 
studied to mimic the transport in Rayleigh-Benard system [11, 12|. The snapshots of the 



velocity field are made of closed lines. The syncronization between the circulation in the 
cells and their global oscillation provides however for a very efficient way of jumping from 



cell to cell. This mechanism, similar to stochastic resonance fl3] , [14]] , makes the curve of 
the effective diffusivity vs the frequency ui of oscillation of the flow very structured. In 
the limit where the molecular diffusion vanishes, anomalous superdiffusion takes place for 
narrow windows of values of uj around the peaks. The anomaly is observed both in the 
second-order moment of particle dispersion and in the singular behavior of the effective 
diffusivity at high Peclet numbers. In Section 4 we study a 2D system consisting of the 
superposition of steady coherent structures and a random incoherent field with a typical 
correlation time r. The resonance occurs now when the sweeping time of the coherent 
structure is close to r. Conclusions are reserved for the final Section 5. 



2 Parallel flow with a large scale velocity field 



The passive scalar field 9(x,t) obeys (see Ref. |15| ) the equation: 



d t 6{x, t) + {v V) 8{x, t) = D A9(x, t). 
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The velocity v(x,t) is here assumed to be incompressible and given by the sum of two 
terms, u and U. The first is periodic both in space (in a cell of size I) and in time 
(the technique can be extended with some modifications to handle the case of random, 
homogeneus and stationary velocity fields). The second is the large-scale component of 
v, varying on a typical scale L. It is assumed that I/L = eCl, where e is the parameter 
controlling the scale separation. 

We are interested in the dynamics of the field 8(x,t) on large scales ~ O(L), com- 
parable to those of U. In the spirit of multiscale methods |J, slow variables X = ex, 
T = e 2 t and r e = et are introduced in addition to the fast variables x and t. The time 
variable r e is needed to account for the advective effects generated by the presence of U. 
The effective equation governing the dynamics at large scales of the field 9 is [IU : 



d t 9 + U-V9 = V a D a pV p 9, (2) 
where the second-order eddy-diffusivity tensor is given by 

D aP (X, T) = 5 al3 D - - \{u a w p ) + {upw a )\ , (3) 

and the symbol (■) denotes the average over the periodicities of the small-scale velocity 
field u. The auxiliary field w(x, t; X,T) has vanishing average over the periodicities and 
satisfies the equation 

d t w + [(u + U) ■ d]w - D d 2 w = -u . (4) 
2.1 A solvable case 

A well-known situation where the effective diffusivities can be calculated analytically is 
the one of parallel flow. Let us consider in detail the case of three-dimensional random 
parallel flow in the presence of a large scale advecting velocity field U(X,T) : 

u(x, t) = (u(y, z, t), 0, 0) ; U(X, T) = (0, U(X, Z, T), 0) . (5) 

The field u(x,t) is random, homogeneus and stationary. The specific dependence on the 
spatial variables makes the fields automatically incompressible. The above decomposition 
is standard in the framework of mesoscale meteorology [[l|J. Its physical meaning is that 
small-scale eddies remain stationary while slow modifications of the large scale component 
occur. This is quite a common feature in geophysics, e.g. in the planetary boundary layer 
]TT| , where the airflow in the thin atmospheric layer near the ground is strongly driven 



by sink/source forcing terms arising from the bottom boundary. 
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The first component of the equation (|j), the only one for which the r.h.s. does not 
vanish, can be easily solved in Fourier space and its solution reads: 



Wi(k, u; X, T) 



-u(k, uj) 



-u(k, u)G(k, u>; X, T) 



iuj + k 2 D + k 
where the advection-diffusion propagator G(k,u; X , T) is defined as: 



(6) 



G(k,u;X,T) 



1 



iuj + k 2 D + iU ■ k 



j+^Do+iU-k]^^ 



It follows from (^|) that the only non-vanishing elements of the eddy-diffusivity tensor 
are diagonal and : 

D XX (X, T)=D + J™ S(t) J E(k)e-^ Do+iU - k ^ dt dk ; D yy = D zz = D , (7) 

We have for simplicity assumed the separability E(k,uS) = E(k)S(u), where E(k) and 
S(uo) are the spectra of the spatial and the temporal part of the correlation function of u. 

The simple analytic expression (|7|) can be used to address the question of interest in this 
paper, i.e., whether it is possible to have a resonance between the oscillation frequencies 
induced by U and the small-scale velocity field. We shall assume for simplicity U to be 
independent of both X and T and the following temporal part of the velocity correlation 
function : 

S(t) = e~ Wr cos(c^) . (8) 
Substituting @ into (0) and performing a simple integration, the latter takes the form: 



Here, the function R(k) is 
R(k) = 



D = D + t E{k)R{k) dk. 



l+k 2 D i 



(9) 



1 _L ( <^L ^ + ( Urk 

~ \ l+k 2 D T J ~ \l+k 2 D T 



Ukr 



,l+k 2 D Tj 1 \l+D k 2 i 

and To is the correlation time of turbulence: 

t = / S(t) dt 



Mki 



(l+D k 2 T) 2 



(10) 



T 



1 + U 2 T 2 



(11) 



Note that R(k) — > 1 when both U — > and D — > 0. In the latter case D = u^tq, 
the well known result valid for velocity fields 5-correlated in time (see, e.g., Ref. |§). It 
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follows from (H) that the combined effect of D and U is to modify the correlation time 
of turbulence, producing an effective correlation time r e ff = r R(k), which depends on 
the wavenumber k. It is precisely such a dependence on k that makes it possible to have 
resonances. Let us indeed consider ([10| ) for a vanishing value of D . It is easy to verify 
that, for ujt > l/y/3, R(k) has a peak for 

1 



jjres = J_ U 2 T 2 _ I + 2uT\/u 2 T 2 + 1. (12) 

kr 

In the limit ujt ^> 1, the previous expression reduces to: U res = Note that R{k) 
crosses the unit value for ujt > 1/ y/3. The regime R(k) > 1 corresponds to constructive 
interference between the turbulent motion and the slow-varying one [|^|. Destructive 
interference takes place, on the contrary, for ujt < l/y/3 where R{k) monotonically de- 
creases from 1 (U — 0) to zero (large U). The effective diffusivity for small-scale flows 
with a single wavenumber ko, attains at the maximum the value 

rye. » (1+^T 2 )(1 + 2^V) 

D = U ° T ° ITI^ ■ (13) 

When ujt ^> 1 (and the velocity U is 'synchronized' with the oscillating velocity field 
(jjres _ ^Ika)), the previous expression reduces to UqTq^-. For U ^> uj/ko, the effective 

diffusivity D goes to zero quadratically, as D = u^tq (j^j) ■ 

The final result is that for a random parallel flow periodic in a box of size L (fco = 2tt/L) 
and ujt ^> 1, the peaks of the effective diffusivity D occur for 

K es = ^ n=l,2,..., (14) 

with values given, for large U, by : 



uj 2 t 2 



~ E(k n)T — . (15) 

The behavior of D as a function of U, obtained evaluating numerically the integral (|9]), is 
shown in Fig. ffl for a Kolmogorov spectrum E{k) cx k~ 5 ^ 3 , and for three different values 
of ujt > l/v^37 Peaks of enhanced diffusion are evident for values of U given by ([14]). 
Note that for ujt = 40, the first five peaks of D (from right to left) are close to the values 
given by the expression (|i~5|). 

The role of a non-vanishing molecular diffusivity D Q is to reduce the height of the 
peaks. This point is shown in Fig. 0, where the behavior of D/(uqTq) is plotted as a 
function of UkoT, for ujt = 40 and for three different values of Dq. As we shall see in 
the following, these features hold also for more general flows than those considered in this 
section. 
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Figure 1: The eddy-diffusivity D/(uqTq) vs the mean flow strength Uk^r. The curve is 
obtained by numerical integration of @ for Dq = 0, E{k) oc k~ 5 ^ 3 , temporal correlation 
function S(t) = e~^^ T cos(tot) and Strouhal number S = MoTo^o ~ 1- The three curves 
correspond to (a): cut = 40; (b): tor = 5; (c): cut = l/y/3. The dashed line joins the 
peaks of D given by the expression (|l5) for n = 1, 2, • • • , 5. 



It is finally worth noting that parallel flow permit to present yet another simple ex- 
ample of resonance. Considering indeed the expression (^]) for vanishing values of U, it is 
easy to verify that for 

ut- 1 u 

Do = n^r~¥ (16) 

peaks of D occur. This corresponds to a synchronization between the characteristic vis- 
cous frequency k 2 D of the random molecular noise and the frequency uo of the oscillatory 
velocity field. 



2.2 Parallel flow with moderate scale separations 

The aim here is to show that the resonant enhanced diffusion previously discussed also 
appears when only moderate scale separation factors are present. As in (jl]), consider 
indeed the velocity field v = u + U with u parallel to the x axis 

u — (ui cos(kxy + ujxt), 0, 0) (17) 

and U parallel to the y axis 

U = (0,u 2 cos(k 2 x + Lu 2 t),0 ) (18) 
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Figure 2: The same as in Fig. [I], but the behavior of D/(uqTq) is now plotted for three 
different values of D and lot = 40. Dotted line: D = 0; full line: -Do/(^o r o) = 5 x 1CT 4 ; 
dashed line: Dq/{uqTq) = 5 x 10~ 3 . Note that the larger is D Q , the broader are the peaks. 



and k\ > k 2 . The situation ki ^> k 2 corresponds to the scale-separated case previously 
considered. Here, we shall assume k\ = 2.5 k 2 (and Ui = u 2 ) and we shall be interested in 
the dynamics at scales much larger than l/k 2 . The two components of the velocity field 
v have moderate scale separation and affect transport on equal footing. The velocity v 
not being a parallel flow, analytical expressions for the effective diffusivities are not 
available anymore. We had then to solve numerically the standard auxiliary equation 
arising from multiscale methods (see, e.g., Ref. [§])■ The equation is solved by using 
a pseudo-spectral method [^] in the basic periodicity cell with a grid mesh of 64 x 64 



points. De-aliasing has been obtained by a proper circular truncation which ensures better 
isotropy of numerical treatment. Time marching has been performed using a leap-frog 
scheme mixed with a predictor-corrector scheme (see Ref. [p2]| ) at regular intervals. 

To enlighten the resonant diffusion enhancement, u\ is fixed and the Dij behavior as 
a function of u> 2 is investigated. Results are reported in Figs. |3] and |[ where Dh/{uqTi) vs 
4u 2 /ui is presented for different values of D . Here, Uq = u\ + u 2 and t\ = 2^ is the 
oscillation period of the small-scale component. Results clearly exhibits resonance peaks 
for uj 2 = ui, u 2 = Ui/2, uj 2 ^ oj sw and uj 2 ~ 2u sw , where the frequency u sw is defined as: 

2tt L 
uj sw = — with T sw = — 

■J-sw 
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Figure 3: The non-dimensional eddy-diffusivity Du/(uqTi) vs the non-dimensional fre- 
quency Aufy/uJii for different values of Dq/{uqT\). The cases (a), (b) and (c) correspond to 
Do/iuln) = 3 x 1(T 3 , D /(uln) = 3 x 1(T 4 and ZV(«o r i) = 3 x 10_f \ respectively 

and 4 u; su ,/co>i = \/2. There are other peaks for u < lu sw , corresponding to sub multiples 
of iO\ and u) sw . 

Note in Fig. | that D 2 2 increases for u 2 — > 0. This is qualitatively understood as follows 
: for small u 2 , the velocity field is essentially a stationary shear flow in the y direction 
superimposed to a rapidly changing small-scale component. The picture is rather close to 
the Taylor mechanism for the longitudinal diffusion in the shear flow, which is well known 



to lead to a 1/D behavior for D 22 p3[ . 



3 Diffusion in time-periodic Rayleigh-Benard convec- 
tion 

Let us consider the two-dimensional velocity field defined by the following stream function : 

ip(x, y, t) = ip sin(x + B smut) sin y . (19) 



This flow [11, |12| is a simple model for transport in time-periodic Rayleigh-Benard 
convection. The stream function (|T9"| ) describes single-mode, two-dimensional convection 
with rigid boundary condition. The even oscillatory instability [lj| is accounted for by the 
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Figure 4: The non dimensional eddy-diffusivity D22I iu^i) vs ^ ne non dimensional fre- 
quency 4co> 2 /^i, for different values of Dq/{uqT\). The cases (a), (b) and (c) correspond to 
Do/(ujjri) = 3 x lO" 3 , D /(«qTi) = 3 x and ZV(«o r i) = 3 x 10" 5 , respectively. 



term B sin cut, representing the lateral oscillation of the roll. In Ref. |fL2j| , a quantitative 
comparison of the behavior in this flow with the experimental data has shown that the 
basic mechanisms of convective transport are well captured by the expression ([19|). 

For a fixed cut, the oscillation frequency of a Lagrangian particle close to the center 
of the periodic cell is of the order of ipQ. The periodicity of the cell is L = 2ir and 
the dimensionless parameter controlling the dynamics is e = ujL 2 The two limiting 
regimes e < 1 and e>l have been investigated in Ref. |2(| to obtain the expressions of the 



diffusion coefficients in the limit of zero molecular diffusivity. Here, we shall concentrate 
on the behavior for e ~ 1. It is indeed precisely in this region that a synchronization 
between the characteristic frequencies O(ipo) and 0(u) can take place. In all the cases 
reported here, a resolution 512 x 512 has been used, which is found adequate for molecular 
diffusivities D /ip > 5 x 1CT 4 . The system evolution has been computed up to times 
500 L 2 /ip , with a time step At = lO^L 2 /^- 

For small enough (or vanishing) numerical diffusivities, the numerical integration of 
the partial differential equation (Q) is clearly subject to instabilities. To evaluate the 
eddy-diffusivities, it is then more convenient to integrate the equation dtx = v using a 
second-order Runge-Kutta scheme and then performing a linear fit of {[x(t) — x(0)] 2 ) vs 
t, to obtain the diffusion coefficient. The averages are made over different realizations 
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oj L 2 / 'i^ 



Figure 5: The turbulent diffusivity D/tpo vs the frequency u)L 2 /tp for different values of 
the molecular diffusivity D /i[) . On the left: D /ipo = 3 x 1CT 3 (dotted line); -Do/^o = 
1 x l(r 3 (dashed line); D /ip = 5 x 10~ 4 (full line). On the right: D /ip = 1 x 1(T 2 (a); 
A)M) = 6 x 10- 3 (b); DoM, = 3 x 10- 3 (c); D /^ = 1 x 1(T 3 (d); £> M> = 5 x 1(T 4 
(e). 

and are performed by uniformly distributing 10 6 particles in the basic periodic cell. The 
system evolution has been computed up to times 10 4 L 2 /ipo- Numerical schemes have been 
tested both for vanishing values of u and for large Peclet number. Asymptotic solutions 
for the effective diffusivity are indeed available in these limits (see Eq. (27) in Ref. p4|). 
Measured diffusion coefficients agree with the theoretical values within errors of less than 
1%. 

The turbulent diffusivity (the x-component) vs the frequency u)L 2 /ijj is shown in Fig. |^ 
for different values of D Q /ijj . A few comments are in order. First, the turbulent diffusivity 
shows maxima which seems similar to those observed for stochastic resonance |T3| , |14| . 
Here, the resonance is between the lateral roll oscillation frequency and the characteristic 
frequencies of the scalar field motion. Second, the effect of Do on the shape of the peaks 
is twofold: similarly to the case of the parallel flow, the smaller is Do, the sharper are 
the peaks. Furthermore, from the expanded view of Fig. [5] (on the right), it appears that 
variations of -Do also cause a shifting of maxima positions. 

For Do — > 0, the peak shown in Fig. moves toward uL 2 /ipo ~ 1-1 and the transport 
becomes superdiffusive at D = 0, i.e. the mean-squared displacement is superlinear: 

([x(t) -x(0)] 2 ) oc t a with a>l . (20) 
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Figure 6: The mean-squared displacement vs the time for the flow (|19|) with D = (in a 
log-log scale). From above to below: uL 2 /ip = 1.075, uL 2 /ip = 1.1 and uL 2 /tp = 1.15. 
Dashed lines have slopes 1, 1.4 and 1, respectively. 
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Figure 7: The diffusion coefficient D as a function of Dq at ujL 2 /ipo = 1.1 (in a log-log 
scale). The slope of the dashed line is —(3 ~ 0.18. 

The range of uoL 2 /ipo values where the process is super diffusive turns out to be very 
narrow: for u>L 2 /ipo = 1.1 the mean-squared displacement behavior is given by (p0|) with 
exponent a ~ 1.4; both for u;L 2 /ip = 1.075 and uL 2 /ip = 1.15 the diffusion process 
is standard (i.e., a = 1). This point is illustrated in Fig. || where the mean-squared 
displacement as a function of the time is shown in a log-log plot for : u>L 2 /ifjQ = 1.075, 
u)L 2 /ip = 1.1 and u)L 2 /ip = 1.15. 

A more careful test of the anomalous diffusion at D Q = consists in computing the 
diffusion coefficient D as function of -Do- As suggested in Ref. |§, in the presence of 
genuine anomalous diffusion, the effective diffusivity must diverge and it is expected that 
D ~ Dq 13 with (3 > 0. The curve D versus Dq at uoL 2 /tpo = 1.1 is shown in Fig. |7|. The 
data are well fitted by a straight line with slope (3 ~ 0.18, confirming the presence of 
anomalous diffusion at -Dq = 0. 



4 Enhanced diffusion in the presence of coherent struc- 
tures 

We shall consider now the case where the two-dimensional large-scale component U con- 
sists of a coherent structure. The small-scale part u is random. Specifically, u = («i, w 2 ), 
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where the u^s are homogeneous, isotropic and stationary Gaussian random processes. 
Their Fourier transforms have zero mean and correlation functions : 



(ui(k, t)u^hf, if)) = (2n) 2 S(k + k')R i:j {k)T(k, t - if) (21) 

where 

Rijik) = (5ijk 2 - kikj)f(k) and T{k,t - f) = e"W. 
The f(k) function is related to the field energy spectrum E(k) as follows: 

The following expressions for E(k) and r(k) have been considered: 

A k 3 k G [0, k T J 
E(k) = { Bk-i kelk^kU (22) 



C e~ k k > k T M 



and 



r(k) = { Fk-i k G [C kU (23) 




' e' k > k\ A 



where the constants A, C and D, G are chosen to ensure the continuity of E(k) and r(k) 
in k^ and kj^, respectively. The expression for r(k) is meant to mimic 3D turbulence, 
neglecting intermittency and correlations between phases. The Ui(k,t) can be easily 
generated by a linear Langevin equation. The coherent large-scale component has the 
cellular structure 

U — U (cosy, cos x) (24) 

characterized by the typical Lagrangian time T c = tt/(\/2U). /,From the isotropy of the 
field it follows Dn = D22 = D. In Fig. [8] the diffusion coefficient D/(v,qTi) is shown as a 
function of T c /ti, where T\ = r(k^ l ) is a characteristic time scale of the random field and 
u 2 = JdkE(k). 

The above results have been obtained performing direct numerical simulations of the 
auxiliary equation (ffl) (now U = and u has the spectrum defined by (E^)) by using a 
pseudo-spectral method over the periodic box 2tt x 2n with a resolution 64 x 64. 

It is well evident the presence of a peak for T c /ti ~ 7 and Do/( M o r i) = 3-9 x 10~ 4 , 
whose value is D/(uqTi) ~ 0.86. The fact that the peak does not appear for T c /t\ of order 
unity is most likely due just to the freedom in the definition of the characteristic times T c 
and T\. On the contrary, it is important to note that the value of the peak is much larger 
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Figure 8: The non dimensional eddy-diffusivity D/(uqT\) vs the non dimensional La- 
grangian time T c /t\ for D /(uqTi) = 3.9 x 10~ 4 and k^/k^ = 4.5. 

than the one occurring when either the random field or the coherent structure are absent. 
The value D/{uqTi) ~ 0.86 has in fact to be compared with 0.07 and 0.32, respectively. 

Unlike the resonance in the Reyleigh-Benard convection, the syncronization mecha- 
nism produces here a broad peak. This is not surprising as in the present case the random 
field has a wide spectrum of caracteristic times. 

5 Conclusions 

Multiscale techniques have been used to present a series of explicit examples of strong 
diffusion enhancement due to resonance effects between different transport mechanisms. 
The effects appear already from the analytical expression of eddy-diffusivity for random 
flows with a large scale mean flow component, are robust to a reduction of the scale sepa- 
ration and they are also observed in convective transport and when coherent structures are 
present on the large scales. For vanishing molecular diffusivities and for specific values of 
the control parameters, the enhanced diffusion might turn into anomalous super diffusion. 
Let us finally stress in which respect the resonant enhanced diffusion effects discussed in 
this paper differ from other known mechanisms to produce strong transport. For Taylor 
longitudinal diffusion in shear flow, open streamlines producing large effective diffusivities 
are for example present in the snapshots of the velocity field itself. For the convective 
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flow in Section 3, snapshots of the velocity field are made at every time of closed cells. 
Open streamlines, and more generally enhanced transport, have a dynamical origin and 
require a temporal syncronization between different transport mechanisms. 
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